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Abstract. The subject of this paper is the study of the asymptotic behaviour 
of the equilibrium configurations of a nonlinear ly elastic thin rod, as the di- 
ameter of the cross-section tends to zero. Convergence results are established 
assuming physical growth conditions for the elastic energy density and suitable 
scalings of the applied loads, that correspond at the limit to different rod mod- 
els: the constrained linear theory, the analogous of von Karman plate theory 
for rods, and the linear theory. 



1. Introduction and Statement of the Main Result 

A classical question in nonlinear elasticity is the derivation of lower dimensional 
models for thin structures (such as plates, shells, or beams) starting from the three- 
dimensional theory. In recent years this problem has been approached by means 
of T- convergence. This method guarantees, roughly speaking, the convergence of 
minimizers of the three-dimensional energy to minimizers of the deduced models. 
In this paper we discuss the convergence of three-dimensional stationary points, 
which are not necessarily minimizers, assuming physical growth conditions on the 
stored-energy density. In particular, we extend the recent results of [16] to the case 
of a three-dimensional thin beam with a cross-section of diameter h and subject 
to an applied normal body force of order h a , a > 2. These scalings correspond at 
the limit to the constrained linear rod theory (2 < a < 3), the analogous of von 
Karman plate theory for rods (a = 3), and the linear rod theory (a > 3). 

We first review the main results of the variational approach. Let £lh = (0, L) X hS 
be the reference configuration of a thin elastic beam, where L > 0, S C M 2 is a 
bounded domain with Lipschitz boundary, and h > is a small parameter. Without 
loss of generality we shall assume that the two-dimensional Lebesgue measure of S 
is equal to 1 and 

/ X2dx2dx3— / Xsdx2dx3— / X2X3 dx^dx^ — 0. (1-1) 
Js Js Js 

Let f h <E L 2 (f2^;IR 3 ) be an external body force applied to the beam. Given a 
deformation v e W 1,2 (fi/,,M 3 ) the total energy per unit cross-section associated to 
v is defined as 

F h (v) = ^[ W{Vv)dx-^( f h -vdx, 
h JO.H h Jn h 

where the stored-energy density W : M 3x3 — > [0, +00] is assumed to satisfy the 
following natural conditions: 
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(HI) W is of class C 1 on M 3 + x3 ; 

(H2) W(F) = +00 if dct F < 0, W{F) -> +00 if dct F -> 0+; 

(H3) W(iJF) = VK(F) for every G 5*0(3), -F 1 G M 3x3 (frame indifference); 

(H4) W = on 50(3); 

(H5) 30 > such that > Odist 2 (F, 50(3)) for every F G M 3x3 ; 

(H6) is of class C 2 in a neighbourhood of 50(3). 

Here M 3 * 3 = {F G M 3x3 : detF > 0} and 50(3) = {R G M 3x3 : i? T i? = 
Id, deti? = 1}. In particular, condition (H2) is related to non-interpenetration of 
matter and ensures local invertibility of C 1 deformations with finite energy. 

The study of the asymptotic behaviour of global minimizers of T h ', as h — > 0, can 
be performed through the analysis of the r-limit of T h (see [7] for an introduction 
to T-convergence) . To do this, it is convenient to rcscale Slh to the domain il = 
(0, L)xS and to rescale deformations according to this change of variables by setting 

y(x) := v(x!,hx2,hx 3 ) 

for every i£fl, Assuming for simplicity that f h (x) = f h (x\), the energy functional 
can be written as 

T h {v)=J h {y)= [ W(V h y)dx- f f h -ydx, 
Jn Jn 



where we have used the notation 



„ (r, d 2 y 3 y\ 

Vm--=\*v — —)■ 

Let now y h be a global minimizer of J h subject to the boundary condition 

y h (0,x 2 ,x 3 ) = (0,hx 2 ,hx 3 ) for every (x 2 ,x 3 ) G 5. (1.2) 

The asymptotic behaviour of y h , as h — >• 0, depends on the scaling of the applied 
load f h in terms of h. More precisely, if f h is of order h a with a > 0, then J h {y h ) = 
O^), where /3 = a for < a < 2 and f3 = 2a - 2 for a > 2, and y' 1 converges in 
a suitable sense to a minimizer of the Tdimit of the rescaled functionals h~P J h , as 
h — > (see [10, 3, 12, 13, 18, 19]). In particular, it has been proved in [13, 19] that, 
if f h is a normal force of the form h a (f 2 e 2 + f 3 e 3 ), with a > 2 and f 2 , f 3 G L 2 (0, L), 
then 

y h ->x 1 e 1 in W 1,2 (Q,R 3 ). 
In other words, minimizers converge to the identity deformation on the mid-fiber of 
the rod. This suggests to introduce the (averaged) tangential and normal displace- 
ments, respectively given by 

{t^tt / {Vi - xi) dx 2 dx 3 if a > 3, 
\ JS , (1.3) 
fe2(a _ 2) (Vi - xi) dx 2 dx 3 if 2 < a < 3, 

v k(^) ■= J^Z2 J s Vk dx 2 dx 3 forfc = 2,3 (1.4) 

for a.e. X\ G (0, L), and the (averaged) twist function, given by 

1 1 f 

w h {x 1 ) := Jj^j^=i J s (x2y 3 -x 3 y^)dx 2 dx 3 (1.5) 

for a.e. £1 G (0, L), where /u(5) := J s (x 2 + x 3 ) dx 2 dx 3 . As h — >• 0, the sequence 
(u^, DjjfJ, w h ) converges strongly in Vl^ 1,2 to a limit (u, «2, W3, w), which is a global 
minimizer of the functional J a given by the r-limit of h~ 2a+2 J h . If a = 3, the 



CONVERGENCE OF EQUILIBRIA OF THIN ELASTIC RODS 



3 



T-liniit J 3 corresponds to the one-dimensional analogous of the von Karman plate 
functional. For a > 3 the functional J a coincides with the linear rod functional, 
while for 2 < a < 3 the limiting energy is still linear but is subject to a nonlinear 
isometric constraint (see Section 2 for the exact definition of the functionals J a ). 

In this paper we focus on the study of the asymptotic behaviour of (possibly non 
minimizing) stationary points of J , as h — > 0. The first convergence results for 
stationary points have been proved in [14, 15, 17]. We also point out the recent re- 
sults [1, 2] concerning the dynamical case. A crucial assumption in all these papers 
is that the stored-energy function W is everywhere differcntiable and its derivative 
satisfies a linear growth condition. Unfortunately, this requirement is incompatible 
with the physical assumption (H2). At the same time, if (H2) is satisfied, the con- 
ventional form of the Euler-Lagrange equations of J h is not well defined and it is 
not even clear to which extent minimizcrs of J h satisfy this condition (we refer to 
[5] and [16] for a more detailed discussion). 

Following [16], we consider an alternative first order stationarity condition, in- 
troduced by Ball in [5]. To this aim we require the following additional assumption: 

(H7) 3k > such that \DW(F)F T \ < k(W(F) + 1) for every F G M 3 + x3 . 

This growth condition is compatible with (H1)-(H6) (see [5]). 

Definition 1.1. We say that a deformation y G W 1 ' 2 ^, R 3 ) is a stationary point 
of J h if it satisfies the boundary condition y(Q 1 x 2 ,x 3 ) = (0, hx 2 , hx 3 ) for every 
(x 2 , x 3 ) G S and the equation 



for every G C^R 3 ,^ 3 ) such that <p(0, hx 2 , hx 3 ) = for all (x 2 ,x 3 ) G S. 

In the previous definition and in the sequel C^(IR 3 ,IR 3 ) denotes the space of C 1 
functions that are bounded in R 3 , with bounded first-order derivatives. 

Assuming (HI) (H7) and using external variations of the form y + e<j)oy 1 one 
can show that every local minimizer y of J" 71 , subject to the boundary condition 
y(0, x 2l x 3 ) — (0, hx 2l hx 3 ) for every (x 2 , x 3 ) G S, is a stationary point of J h in the 
sense of Definition 1.1 ([5, Theorem 2.4]). Moreover, when minimizcrs are invertiblc, 
condition (1.6) corresponds to the equilibrium equation for the Cauchy stress tensor. 

In [16] it has been proved that stationary points in the sense of Definition 1.1 
converge to stationary points of the T-limit J a in the case of a thin plate and for 
the scaling a > 3 (corresponding to von Karman and to linear plate theory). In 
this paper we extend this result to the range of scalings a > 2 in the case of a thin 
beam. Our main result is the following. 

Theorem 1.2. Assume that W satisfies (H1)-(H7). Let f 2 , f 3 G L 2 (0,L), and 
a > 2. For every h > let y h be a stationary point of J h (according to Defini- 
tion 1.1) with f h :— h a (f 2 e 2 + f 3 e 3 ). Assume there exists C > such that 




(1.6) 




(1.7) 



for every h > 0. Then, 



y h ^ Xl ei inW h2 (n,R 3 ). 



(1.8) 
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Moreover, let u , v , and w h be the scaled displacements and twist function, intro- 
duced in (1.3)-(1.5). Then, up to subsequences, we have 

u h ^u tnW 1,2 (0,L), 

v%^v k inW 1>2 (0,L) fork = 2,3, 

w h ^w inW 1>2 (0,L), 

where (u,v 2 ,v 3 ,w) £ W 1,2 (0, L)xW 2 < 2 (0, L)xW 2 < 2 (0, L)xW 1,2 (0, L) is a station- 
ary point of J a . 

The proof of Theorem 1.2 is closely related to [16] and uses as key tool the rigidity 
estimate proved in [9]. The main new idea with respect to [16] is the construction 
of a sequence of suitable "approximate inverse functions" of the deformations y h 
(see Lemma 2.7), which allows us to extend the results of [16] to the range of 
scalings a 6 (2, 3). This construction is based on a careful study of the asymptotic 
development of the deformations y h in terms of approximate displacements and 
uses in a crucial way the fact that the limit space dimension is equal to one. 



2. Preliminary Results 

In this section we recall the expression of the F-limits J a identified in [13] and 
[19] and we prove some preliminary results. 

We start by introducing some notation. Let Q 3 : M 3x3 — > [0, +00) be the qua- 
dratic form of linearized elasticity: 



Q 3 {F) := D 2 W(Id)F : F for every F e 



(r3x3 



We will denote by C the associated linear map on M 3x3 given by C := D 2 W{Id). 
Let 

E:= min Q 3 (ei\a\b), (2.1) 

and let Q\ be the quadratic form defined on the space M 3 ^ of skew-symmetric 
matrices given by 

Qi(F):= min / Q 3 ( x 2 Fe 2 + x 3 Fe 3 d 2 [3 8 3 p) dx 2 dx 3 (2.2) 

I3GW 1 ' 2 (S,R 3 ) J s V / 

for every F € M 3 ^. It is easy to deduce from the assumptions (H1)-(H6) that E 
is a positive constant and Q\ is a positive definite quadratic form. 
The functional J a are defined on the space 

H := W 1 ' 2 (0,L)xW 2 ' 2 {0,L)xW 2 ' 2 {0,L)xW 1 ' 2 {0,L) 

and are finite on the class A a , which can be described as follows: 

A a := {(u,v 2 ,v 3 ,w) e H : u' + \[{v 2 ) 2 + (v 3 ) 2 } = in (0, L) and 

u(0) = v k {0) = v' k {0) = w(0) = for k = 2,3} 

for 2 < a < 3, and 

A a := {(u,v 2 ,v 3 ,w) G H : u(0) = v k (0) = v' k (0) = w(0) = for k = 2,3} 
for a > 3. 

For 2 < a < 3 the functional J a is given by 

J a {u,v 2 ,v 3 ,w) = ^ / Qi(A')dxi- [ (f 2 v 2 + f 3 v 3 ) dxi (2.3) 
* Ja Jo 



CONVERGENCE OF EQUILIBRIA OF THIN ELASTIC RODS 



5 



for every (u, V2, v 3 , w) G A a , J a {u,V2,v 3 ,w) = +00 otherwise in H. In (2.3) the 
function A G W^ 2 ((0, L), M 3x3 ) is defined by 

/ -v' 2 (xi) -v' 3 (xi) 

A(xi) := v' 2 { Xl ) -w{xi) \ (2.4 

\ v' 3 (xi) w(xi) 

for a.e. x\ G (0, L). 

For a = 3 the r-limit is given by 

1 

Mu,V2,v 3 ,w) = - J E(«' + |[(^) 2 + (^) 2 ]) dan 



+ 



^jf Q^dxt-J (f2V2 + f 3 v 3 )d Xl (2.5) 



for every (w,v 2 ,V3,u>) <E J 3 (u,V2,v 3 ,w) = +00 otherwise in if. 
Finally, for a > 3 the r-limit is given by 

1 /" L 1 /" L 

Ja(«,«2,u 3 ,tu) = g J ®{u') 2 dx 1 + -J Q 1 {A')dx 1 

- f (f2V2 + f 3 v 3 ) <tei (2-6) 
Jo 

for every (w, V2, V3, u>) 6 A a , J a (u,V2,v 3 ,w) = +00 otherwise in H. 

We can now compute the Euler-Lagrange equations for the functionals J a intro- 
duced above. We first recall the following lemma. 

Lemma 2.1. Let F G M 3 ^ and let Q F : W 1 ' 2 ^, K 3 ) -> [0, +00) be the functional 

Qf{P):= J Q 3 (x 2 Fe2 + x 3 Fe 3 d 2 [3 d 3 [3) dx 2 dx 3 

for every (3 G W 1 ' 2 (S,R 3 ). Then Q F is convex and has a unique minimizer in the 
class 

B := |/3 G W 1 ' 2 (S,M 3 ) : J /3 dx 2 dx 3 = J d 2 /3 dx 2 dx 3 = J d 3 /3dx 2 dx 3 = o|. 

Furthermore, a function (3 G B is the minimizer of Qf if and only if the map 
E:S^ M 3x3 defined by 



d 2 (3 



0: 



E := C[x 2 Fe2 + x 3 Fe 3 
satisfies in a weak sense the following problem: 

j div X2 ^ 3 (Ee 2 \Ee 3 ) = in S, 

\{Ee2\Ee 3 )u ds = on dS, 

where vqs is the unit normal to dS . Finally, the minimizer depends linearly on F . 

Proof. See [14, Lemma 2.1] and [12, Remark 3.4]. □ 

We shall use the following notation: for each F G L 1 (J7,M 3x3 ) we define the 
zeroth order moment of F as the function F : (0, L) — > M 3x3 given by 



F{x\) := / F(x)dx2dx 3 
Js 



for a.e. X\ G (0, L). We also introduce the first order moments of F as the functions 
F, F : (0, L) -> M 3x3 given by 

F(x\) := / X2F(x) dx2dx 3 , F(xi) = / x 3 F(x)dx2dx 3 
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for a.e. x\ G (0, L). 

The following proposition follows now from straightforward computations. 

Proposition 2.2. Let (u, v 2 ,v 3 ,w) G A a . For a.e. x\ G (0, L) let /3(x\, •, •) G B be 
the minimizer of Qa'( Xi ), where A' is the derivative of the function A introduced in 
(2.4). Let also E : ft -> M 3x3 be defined by 



E := c(x 2 A'e 2 + x 3 A'e 3 | d 2 /3 | d 3 /?) 



and /ei -E 1 and 6e zis first order moments. Then 

(1) (u,v 2 ,V3,w) is a stationary point of J 3 if and only if the following equations 
are satisfied: 

«' + \[W 2 ? + (v' 3 f] = in(0,L), (2.8) 



(2.9) 
(2.10) 
(2.11) 



te+/2=0 m(0,L), 
\fii(L) = ^ 11 (L)=0, 
f^i+/3 = m(0,L), 
\Sn(L) = ^ 1 (L)=0, 
r^ 2 = ^3 m(0,L), 
\e 12 (L) = E 13 (L); 

(2) if a > 3, (u,v 2 ,v 3 ,w) is a stationary point of J a if and only if 

u' = m(0,L) (2.12) 

and (2.9)- (2.11) are satisfied; 

(3) if 2 < a < 3, (u,v 2 ,v 3 ,w) is a stationary point of J a if and only if (2.9)- (2.11) 
are satisfied. 

Remark 2.3. If (u,v 2 ,v 3 ,w) G A a and 2 < a < 3, then u is uniquely determined 
in terms of v 2 and v 3 . Indeed by the constraint 

u / + ^) a + K) a =0a .e. i n(0,L) 
and the boundary condition u(0) = 0, we have 

r Xi (i/ ~) 2 -+- (v' ~) 2 

u(xi) = -/ y 2 ' ^ { - for a.e. x x in (0,L). (2.13) 

For a > 3 the same conclusion holds when (u, v 2 ,v 3 ,w) G -4 a is a stationary point 
of J a . Indeed, if a = 3, (2.8) yields (2.13), while, if a > 3, (2.12) gives 

u = a.e. in (0, L). 

Using the previous observations and the strict convexity of Q\, it is easy to show 
that for every a > 2, J a has a unique stationary point that is a minimum point. 

Remark 2.4. For what concerns the three-dimensional functionals J h , under ad- 
ditional hypotheses on W (such as polyconvexity, see [4]) it is possible to show 
existence of global minimizers, and therefore of stationary points. Furthermore, 
they automatically satisfy the energy estimate (1.7) (see [10, proof of Theorem 2]). 
For general W the existence of stationary points (according to Definition 1.6 or 
to the classical formulation) is a subtle issue. We refer to [5, Section 2.7] for a 
discussion of results in this direction. 
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From now on we shall work with sequences of deformations y h e W 1,2 (fl, R 3 ), 
satisfying the boundary condition (1.2) and the uniform energy estimate (1.7) with 
a > 2. This bound, combined with the coercivity condition (H5), provides us with a 
control on the distance of Vhy h from SO(3). This fact, together with the geometric 
rigidity estimate by Frieseckc, James and Miiller [9, Theorem 3.1], allows us to 
construct an approximating sequence of rotations (R h ), whose L 2 -distance from 
\7hU h is of the same order in terms of h of the L 2 -norm of dist(V/ l y' 1 , SO (3)). More 
precisely, the following result holds true. 

Theorem 2.5. Assume that W : M 3x3 — > [0, +oo] is continuous and satisfies 
(H3)-(H6). Let a > 2 and let (y h ) be a sequence in W 1 ' 2 {9., 1 M?) satisfying (1.2) 
and (1.7) for every h > 0. Then there exists a sequence (R h ) in C°°((0, L), M 3x3 ) 
such that 



R h (x 1 ) E 50(3) for every Xl e (0,L), (2.14) 

\\V h y h - R h \\ L 2 <Ch a -\ (2.15) 

\\(R h y\\L* <Ch- 2 , (2.16) 

\\R h -Id\\ L ~ <Ch a -\ (2.17) 

We omit the proof as it follows closely the proof of [15, Proposition 4.1]. Owing 
to the previous approximation result, one can deduce the following compactness 
properties. 

Theorem 2.6. Under the assumptions of Theorem 2.5, let u h , v\, v^, w h be the 
scaled displacements and twist function introduced in (1.3)-(1.5). Then 

y' l ^xiei strongly in W 1,2 (0,R 3 ) (2.18) 

and there exists (u,V2,Vs,w) £ A a such that, up to subsequences, we have 

u h ^u strongly in W h2 (0,L) if2<a<3, (2.19) 

u h u weakly in W 1 ' 2 ^, L) if a > 3, (2.20) 

v'^v k strongly in W h2 (0,L), fc = 2,3, (2.21) 

w h -± w weakly in W 1 ' 2 ^, L). (2.22) 

Moreover, let A G T4 7l,2 ((0, L), M 3x3 ) be the function defined in (2.4). Then, if 
R h is the approximating sequence of rotations given by Theorem 2. 5, the following 
convergence properties hold true: 

^hV h -Id yA ^ r2 . 3x3 



h a ~ 2 



-> A strongly in L (Q,M ), (2.23) 



uh _ T J 

A h := n A weakly in W 1>2 ((0, L), M 3x3 ), (2.24) 

sym(R h -Id) A 2 f , . T ^ 

l 2[a _ 2) ' -» -y uniformly m (0, L). (2.25) 

For the proof we refer to [19, Theorem 3.3]. 



We conclude this section by proving a lemma, which will be crucial to extend 
the convergence of equilibria result to the scalings a e (2, 3). 
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Lemma 2.7. Under the assumptions of Theorem 2.5, there exist two sequences 
(£fc)> & = 2, 3, such that for every h > 



$eCt(R), ££(o) = o, 



Vk 



1 



- -^£k ° Vi -> ^fc strongly in L 2 (n), 

I^II^ + II(^)'IIl~<c^- 2 . 



(2.26) 

(2.27) 
(2.28) 



Remark 2.8. The sequences of the previous lemma can be interpreted as 
follows: the functions defined by 

u [x) ~ r 1 ' h h ' h h ) 

represent a sort of "approximate inverse functions" of the deformations y h , in the 
sense that the compositions uj h oy h converge to the identity strongly in L 2 (f2,IR 3 ) 
by (2.18) and (2.27). 

Proof of Lemma 2.7. In order to construct the functions we first study the 
asymptotic behaviour of the sequences k = 2, 3. By Poincare inequality we 

obtain the estimate 



^ - x k - J - a*) da^da*^ < c( 



d k y h k 



h 



- 1 



L 2 



+ 



h 



L- 



)■ 



where k,j G {2,3}, fc ^ j. Therefore, by (1.4) and (2.23) we have 

< ch a - 2 . 



$ - n - h a - 3 v h 



L 2 (Q) 



(2.29) 



In particular, for a > 3 it follows that y£ — > xj, strongly in L 2 , so that, if a > 3, 
we can simply take = for A; = 2, 3 and every ft > 0. If 2 < a < 3, we need to 
construct a suitable approximation of v\. Let (i? ft ) be the approximating sequence 
of rotations associated with (y h ) (see Theorems 2.5 and 2.6). By (2.17) and (2.25) 
we deduce the following estimates: 

< Ch a - 2 forfc = 2,3, 



111 Too KCh 2 ^. 



L ll ~ -Hli 

Let r£,ri G C(K) be continuous extensions of the functions and i?^ — 1 to 
such that for every h > 

supp r£, supp r\ C (-1,L + 1) 



in (0, L) for fc = 2, 3, 
<-lin (0,L), 



\\rkh°° < Ch a - Z forfc = 2,3, 
llrfHio. <Ch 2{a - 2 'K 
We introduce the functions v^, u£ G C^(R) defined by 



v*( Xl ) := r 
Jo 



4{s)ds, 



r?(s)ds. 



(2.30) 
(2.31) 
(2.32) 
(2.33) 
(2.34) 

(2.35) 
(2.36) 



Using the boundary condition (1.2), the Poincare inequality, (2.15), and (2.23), we 
obtain 

1 



h Xk h Vk 



L 2 



< Ch° 



(2.37) 
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and analogously, 

\\v$-xi-t\\t*<Ch a - x . (2.38) 
This last inequality, together with (2.34), implies that 

hi -xi\\ L 2 < Ch 2{a - 2) for a < 3. (2.39) 

We are now in a position to construct the maps £ k when a < 3. If a = 3, we 
define = v%. Properties (2.26) and (2.28) follow immediately. To verify (2.27) it 
is enough to remark that by (2.33) and (2.37) we have 

y£ j} h k oy>? . 1 



l2 < Ch+-\\v^oy^-v h k \\ L2 

< Ch + ^\\(v h k y\\ L ^\\y h 1 -x 1 \\ L2 <Ch. 
If 2 < a < 3, we first fix n E N such that 

" >2+ 2^T3 (2 ' 40) 
and we introduce a sequence of maps n — 1, . . . , n , recursively defined as 

(n(xi) = xx -v^o( n+1 ( Xl ) for n = l,...,n - 1. 
For k — 2, 3 and every ft, > we define 

e k :=v h k oQl (2.42) 

Since C^ (0) = 0, we have by induction that C«(0) = f° r cacri n = 1)2,..., no, so 
that £^(0) = 0. From the regularity of and v k it follows that (2.26) is satisfied. 
By (2.33) we deduce 

< <(L + 2)\\r h k \\ L ~ < Ch-\ (2.43) 

To estimate ||z,°°, we first deduce a recursive bound for ||(Cn)'IU°°- If h is 

small enough, we have 

ll(«l fc )'l|L~ < I- 

By (2.41) the following inequalities hold true: 

\\(C y\\L~<i+\\(vmL~<2, (2.44) 

IKC^'lUoo <l + ||(C5 + i)'IU~ forn = l,...,no-l, (2.45) 

||(Ci)'IU~ < 1 + no- (2.46) 
Now by (2.46) and (2.33) we have 

< ||(^)'|k~||(Cf)'IUoo < (l + n )\\r h k \\ L ~ < Ch<*-\ (2.47) 

Combining (2.43) and (2.47) we obtain (2.28). To conclude the proof it remains to 
verify (2.27). By (2.34), (2.38), and (2.39) we have 



HCi°2/?-Zi|b = hi-v^oy^-x^ 



< lli/f-ef-^ii^ + iief-sfo^ii^ 

< Ch"- 1 + h 2{a - 2) \\r^\\ L oo < Ch*- 1 + Ch A{a - 2) . (2.48) 
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Arguing analogously for Cn -i an( i using (2.48), we obtain 

< cho-' + wivtyu-wCoyt-x^ 

< Ch"- 1 + ch^-^ih"- 1 + h 4 ( a - 2 ')) 

< Ch a - X +Ch 6(a - 2) . (2.49) 

By induction we deduce 

\\(n°yi -Xl\\ L 2 < Ch"- 1 +C -^2(„o-™+2)(a-2)^ y ^ 

In particular, we have 

HCi °Vi -xih* < Ch - 1 + Ch^ na+1 ^ a - 2 l (2.51) 
We can now prove (2.27). By (2.42), (2.33) and (2.51) we obtain 

l\\ti o y\ v h k \\ L , = ±\\v h k o d fc o y\ v h k \\ L , 

<^ll(^)'IU~IICN^-^ilU 2 

< i||^||L=o(C/ l a - 1 + C/i 2( " 0+1)(a - 2) ) 

< Ch a - z (Ch a - x + c/ l 2(«o+i)(«-2)^ 

< £.^min{2a-4 : (2n +3)a-(4ri + 7)} 

where the last term converges to zero because of (2.40). Combining this with (2.37), 
we deduce (2.27). □ 

3. Proof of the Main Result 

This section is entirely devoted to the proof of Theorem 1.2. The proof strategy is 
similar to [16] . The major difference is in the analysis of the asymptotic behaviour of 
the first-order stress moments (Steps 6 and 7), where the approximating sequences 
constructed in Lemma 2.7 are needed to define suitable test functions in the scalings 
2 < a < 3. 

Proof of Theorem 1.2. Let (y h ) be a sequence of deformations in W 1 ' 2 ^, R 3 ) satis- 
fying the energy bound (1.7), the boundary condition (1.2), and the Euler-Lagrange 
equations 

DW(V h y h )(V h y h ) T :[(V<f>)oy h ]dx = f h a [f 2 (h oy h ) + f 3 (<j> 3 oy h )] dx (3.1) 

for every (f> G C£(M 3 ,R 3 ) such that 4>(0,hx 2l hx 3 ) = for all (#2,2:3) G 5". 

Convergence of the sequences (y h ), {u h ), {v k ), and (w h ) follows from Theo- 
rem 2.6, together with the fact that (u, V2,v 3 ,w) G A a . To conclude the proof we 
need to show that (u,V2,Vs,w) is a stationary point of J a . 

The proof is split into seven steps. 

Step 1. Decomposition of the deformation gradients in rotation and strain 

Let (R h ) be the approximating sequence of rotations constructed in Theorem 2.5 
and let A G W ll2 ((0, L), M 3x3 ) be the function defined in (2.4). We introduce the 
strain G h : O -> M 3x3 as 

V h y h = R h (Id + h a ^G h ). (3.2) 

By (2.15) the sequence (G h ) is bounded in L 2 (ft,M 3x3 ), so that there exists G e 
L 2 (ft,M 3x3 ) such that G h G weakly in L 2 {Vt, M 3x3 ). Moreover, by Lemma 3.1 
(see the end of this section) the symmetric part of G can be characterized as follows: 
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there exists [3 £ L 2 (0, R 3 ), with zero average on S and dk(3 £ L 2 (Q, R 3 ) for k = 2, 3, 
such that, if we set 

M(/3) := (x 2 A' e2 + z 3 A' e3 | 8 2 [3 | 9 3 /?) , 

we have 

'symM(/3) + (u' + \[{v' 2 ) 2 + (« 3 ) 2 ])ei ® ei if a = 3, 
symG = < symM(^) + u'd <g> e x if a > 3, (3.3) 

symM(/3) + #ei <g> ei if 2 < a < 3 

for some g £ L 2 (0, L). In particular, by the normalization hypotheses (1.1) on S we 
deduce 

\u' + \[{v' 2 f + {v' 3 f] for a = 3, 
Gn = < u' for a > 3, (3.4) 

[g for 2 < a < 3. 

Step 2. Stress tensor estimate 

We define the stress E h : il ^ M 3x3 as 

£' 1 = -J-^DWild + ft a - 1 G h )(Jd + h a - 1 G h ) T . (3.5) 

From the frame indifference of VF it follows that 

£W(.F)F T = F(LW(F)) T for every F £ M 3 ^ 3 . 

This implies that E h is symmetric for every h > 0. Moreover, the following pointwise 
estimate holds: 

|*|< C ( Sg±£^ + |G>|). (3.6) 

Indeed, let <5 be the width of the neighbourhood of 50(3) where 14 7 is of class 
G 2 . Suppose first that /i Q_1 |G' l [ < |. Then, a first order Taylor expansion of DW 
around the identity, together with (H4) and (H5), yields 

DW(Id + h a - 1 G h ) = h a - 1 D 2 W{M h )G h 

for some 

M h e M 3x3 sat i s f y i n g (m 71 - Id\ < §. Since D 2 W is bounded on the set 
{F £ M 3x3 : dist(F,50(3)) < §}, we deduce 

|W(H + r 1 G'')| < Ch a - l \G h \. 

Therefore, by (3.5) we obtain 

\E h \ < C\G h \ + Ch a - 1 \G h \ 2 < C(l + S)\G h \. 

If instead h a ~ 1 \G h \ > f, we first observe that W{V h y h ) is finite a.e. in D, by (1.7). 
By (H2) and by frame indifference we deduce 

det(V h y h ) = det(Id + h a - l G h ) > a.e. in ft. 

Therefore, we can use (H7), which yields 



\ E »\ < ——k(w(id + hr-W) + 1) < k w ( Id + ha lGh ) + ™\g*\ 

This completes the proof of (3.6). 

Step 3. Convergence properties of the scaled stress 

Arguing as in [16], some convergence properties of the stresses E h can be deduced 
from (3.6). Indeed, using (1.7) and the fact that the G h arc bounded in L 2 (ft, M 3x3 ), 
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we obtain from (3.6) that for each measurable set A the following estimate holds 
true: 



L 



\E h \dx < Ch"- 1 + C|A|*, (3.7) 

where |A| denotes the Lebesgue measure of A. Let now 

B h := {x e ft : h^-'^ix)] < 1}, (3.8) 

where 7 G (0, a — 2), and let \h be the characteristic function of Bh- By (3.7) and 
by Chebyshev inequality we have 

\E h \dx < Ch 01 - 1 -', (3.9) 

n\B h 
so that 

(1 - Xh)E h -> strongly in L 1 (ft, M 3x3 ). (3.10) 

Moreover, one can show that the remainder in the first order Taylor expansion of 
DW(Id + h a ~ 1 G h ) around the identity is uniformly controlled on the sets B h , so 
that 

Xh E h ^CG=:E ini 2 (ft,M 3x3 ) (3.11) 
(see Step 3 in the proof of [16, Theorem 3.1] for details). 

Step 4. Consequences of the Euler- Lagrange equations 
By the frame indifference of W and by (3.2) we have 

DW(V h y h ){V h y h ) T = h a - 1 R h E h (R h ) T . 

Therefore, the Euler-Lagrange equations (3.1) can be written as 

/ R h E h (R h ) T :[(V4>)o y h ]dx = h [ [h{4>2oy h ) + h{fooy h )]dx (3.12) 

for every <f) E C\ (M 3 , K 3 ) satisfying the boundary condition <p(0,hx2,hxs) = for 
all (x2, X3) e S. 

Let now be a function in C^(R 3 ,M 3 ) such that (j)(0,X2,x 3 ) = for every 
(x2,xs) e S. For each h > we define 

4>\x) := W{xu f - J^i), f - , 

where £2; £3 are the functions constructed in Lemma 2.7. By (2.26) the maps (j> h 
are admissible test functions in (3.12). 

To simplify computations we introduce the following notation: 

h. { 'h V2 ch„ n ,h V3 ^ f/i „„,/A fo 1 q\ 

From (1.8) and (2.27) it follows that 

z h ^x inL 2 (ft,M 3 ). (3.14) 
Choosing <f> h as test function in (3.12) we obtain 



r r 3 i 

/ R h E h (R h ) T e 1 -[hd 1 cj>oz h -Y,(dk<i>°z h )(($y°y'?) </■< 

"' n k=2 
3 



■ f Y,R h E h {R h ) T e k -{d k ^oz h )dx+ f h 2 [f 2 (4> 2 oz h ) + f 3 (cf> 3 oz h )]dx = Q. 

(3.15) 
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By (3.7) and (2.28) we have 



r r 3 

/ R h E h (R h ) T e 1 - ft^o* fc -V(0 fc 0o 

J ° 1 k=2 



oy h A 



dx 



< 



6 

c\\E h \\ L , (\\hd l( t>\\ L oo + ]T IK^'lUo.) < c(h + h a -% 



k=2 



therefore the first integral in (3.15) converges to zero. Analogously, since f k G 
L 2 (0,L) for fc = 2,3 and <j> k G C^R), the last integral in (3.15) tends to zero. We 
deduce that the second integral in (3.15) must also converge to zero. On the other 
hand, this term can be written as 



J2 R h E h (R h ) T e k ■ (d k 4>o z h ) dx 

'■ k=2 

3 

= / ^Xh^^^-^O 
^ n k=2 

3 

+ / J2(l-Xh)R h E h (R h ) T e k -(d k <l>oz h )dx. 



' fc=2 



By (3.14) and by the dominated convergence theorem we have 

d k <t>oz h -> d k <t> inL 2 (n). 
Thus, by (3.11) and by the fact that R h ->■ Id in L°°(0,L) we deduce 

/ y jX hR h E h (R h ) T e k -(d k <j)oz h )dx^ [ Y^Ee k -d k <t>dx, 



(3.16) 



(3.17) 



' k=2 



k=2 



while by (3.10) we have that the last term in (3.16) tends to zero. We conclude that 



3 

/ ^ Ee k ■ d k <fi dx = 

"' n k=2 



(3.18) 



for every <j) G C£(R 3 ,]R 3 ) such that </>(0,x 2 ,x 3 ) = for all (x 2 ,x 3 ) G S. Therefore, 
the following equations hold true a.e. in (0, L): 

(div X2tX3 (Ee 2 \Ee 3 ) = in S, 
\{Ee 2 \Eez)v as = on 95, 

where i/gg is the unit normal to dS. Moreover, for a.e. x\ G (0, L) 

I Ee k dx 2 dx 3 = for k = 2, 3. (3.20) 
Js 

We conclude that Ee 2 = Eez = a.e. in (0, L) and since E is symmetric, 

£7 = ® ei. 



Step 5. Zeroth moment of the Euler- Lagrange equations 

We now identify the zeroth order moment of the limit stress E. Let ip be a function 
in C^QR) such that ^(0) = 0. We define 

4>{x) = V>(zi)ei. 
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Using <f> as a test function in the Euler-Lagrange equations (3.12) we have 

(R h E h (R h ) T ) n (^' o y'l) dx = 0. (3.21) 



To pass to the limit in the previous equation, we split f2 into the sets Bh and fl\Bh, 
so that we obtain 

/ X h(R h E h (R h f) 11 (ib'°yt)dx+ f (l- X h)(R h E h (R h ) T )u(i>'oyi?)dx = 0. 
Jn Jn 

(3.22) 

By (1.8) and by the continuity of ip' it follows that tp' o y\ converges to ip' in L 2 (fl). 
Therefore, by (3.10) and (3.11) we can pass to the limit in (3.22) and we deduce 

/ E 11 ip'dx 1 = [ E 11 ip'dx = 
Jo Jn 

for every ip <E C^(R) such that ip(0) = 0. This implies that E = £ n ei ® ei = a.e. 
in(0,L). 

Since by frame indifference CH = for every skew-symmetric H G M 3x3 , we 
obtain that £symG = CG = E = 0. The invertibility of C on the space of sym- 
metric matrices yields that symG = 0. Together with (3.4), this implies (2.8) for 
a = 3, (2.12) for a > 3, and g = a.e. in (0, L) for 2 < a < 3. Moreover, by (3.3) 
we deduce that 



sym(0 



d 2 pJ dx 2 dx 3 



d 3 pJ dx 2 dx 3 ) = 0, 



so that, if we introduce j3 : ft — > R 3 defined by 

P : = (Pi,@2 -x 3 J d 3 (i 2 dx 2 dx 3l j3 3 -x 2 J d 2 p 3 dx 2 dx 3 ), 

we have that f3(x\, -, ■) € B for a.e. x\ e (0, L) and 

symG = sym^2^4'e2 + x 3 A'e 3 d 2 (3 d 3 pfj . 



In particular, we have the following characterization of E: 



E = £symG = C[x 2 A'e 2 + x 3 A'e 3 



ds/3). 



d 2 (i 



Since E satisfies (3.19), we deduce from Lemma 2.1 that j3 is a minimizer of the 
functional 

Q A .{p) = J Q 3 (x 2 A'e 2 + x 3 A'e 3 | d 2 (i | 9 3 /?) dx 2 dx 3 . 
In other words, p3 satisfies 

Qi(A')= J Q 3 (x 2 A'e 2 +x s A'e 3 



d 2 (3 



d 3 (3 ) dx 2 dx 3 



(3.23) 



for all a > 2. 



Step 6. First-order moments of the Euler-Lagrange equations 

In this step we prove that the limiting Euler-Lagrange equations (2.9) and (2.10) 
are satisfied. Let (p 2 ,<p 3 be two functions in C£(R) with ^2(0) = <p 3 (0) = 0. We 
define 

<P2{xi) f 3 (xi)' 



A*)=(o,^,^) 



CONVERGENCE OF EQUILIBRIA OF THIN ELASTIC RODS 



15 



and we use cj> h as test function in (3.12). By (1.8) the force term can be treated as 
follows: 

lim f h[f 2 {4> h 2 o y h ) + h{4> h z o y h )\ dx = lim f [f 2 (<p 2 o j,*) + f 3 (<p 3 o y*)] dx 
Jn Jn 



-I 

Jo 



(/2P2 + hvsjdx-L. 



(3.24) 



Therefore, we have 



lim 

ft— s-0 



(R h E h (R h ) T ) 2l ^^- + (R h E>\R h f) 31 ^±- 



f 

Jo 



h 



(/2V2 + hvsjdxi. 



dx 



(3.25) 



We shall characterize the limit on the left-handside of (3.25) in terms of the first- 
order moments of the stress E. To this aim, we go back to the Euler-Lagrange 
equations (3.12) and we construct some ad-hoc test functions with a linear be- 
haviour in the variables x 2 , x 3 . Let (uih) be a sequence of positive numbers such 
that 



htOh — > +00, 

h^-^LOh -> 0, 



(3.26) 
(3.27) 



where 7 e (0, a — 2) is the same exponent introduced in (3.8). For each h > we 
consider a function 6 h G C\ (R) which coincides with the identity in a large enough 
neighbourhood of the origin, that is, 



9 h (t) = t for \t\ < u h 
and, in addition, satisfies the following properties: 

\9 h {t)\ <\t\ VieR, 



d6 h 



dt 



< 2to h , 



< 2. 



(3.28) 



(3.29) 
(3.30) 

(3.31) 



Let i) be a function in C 1 (R) with compact support and such that 77(0) = 0, and 
let , k = 2, 3, be the functions constructed in Lemma 2.7. We consider the map 



<f> h (x) = 9 h (^-^( Xl ))r 1 (x 1 )e 1 . 

Choosing <ft h as test function in (3.12) and using the notation introduced in (3.13), 
we obtain 

/ (R^^R^n^oz^'oyDdx 
Jn 

- 1 {RhEh f h)T)l1 ° 4) mr o y^Kv o 1A) dx 

+ JjR^E h (R h ) T ) 13 ^{^ozl)dx = 0. (3.32) 
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The first integral in (3.32) can be decomposed into the sum of two terms 

/ {R h E h {R h ) T ) ll {e h oz^){r ] 'oy t l)dx 
Jn 

= f X h[(R h E h (R h ) T )ll(O h o zl)( V ' oy1)]dx 
Jn 

+ I (l- X h)\(R h E h (R h ) T )u(6 h ozl)(r 1 'oy^]dx. (3.33) 
Jn 1 J 

By (1.8), (3.14), (3.29), and by the dominated convergence theorem we deduce that 

(e^zDiv'oy'D^xW in L 2 (n). 
Therefore, by (3.11) we have 



Jim / X h{R h E h {R h ) T ) 11 {i 1 'oy h l ){e h oz%)dx= / x 3 E llV 'dx 



Eur]' dxi. 



The second term in (3.33) can be estimated using (3.9), as follows: 

/ (l- Xh )\(R h E h (R h ) T )u(v'oy^(O h ozl)\dx < 2u h \\rf\\ L ~ {0iL) / \E h \ 
Jn Jn\B h 

< Ch a - 1 -«Lo h , 

and this latter is infinitesimal owing to (3.27). We conclude that 

/ (R h E h (R h ) T ) ll (O h oz*)(r 1 l oy h 1 )dx^ t E n v' d Xl . (3.34) 
Jn Jo 

As for the second integral in (3.32), we consider the following decomposition: 
jf (R h E h (R h f)ul o z*) [«*)' o o y?) dx 

= JjR h E h (R h f) n [(^ozl) -l\lmyoy^}( V oy1)dx 

+ Ji&E^&f^^yoyfitooyftdx. (3.35) 

To study the first term in (3.35) we introduce the sets 

D h = {xeSl: \z$(x)\ >cj h }. (3.36) 
Since (z^) is uniformly bounded in L 2 (tt), by Chebyshev inequality we deduce that 



\D h \ < Cu- h \ 



(3.37) 



Thus, by (2.28) and (3.7) we have 
\( deh „„h 



J^(R h E h (R h f)n [(^ o **) l] oy^oyl) dx 

< j D | ( i?^(i?^) 11 [(^oz^)-l]I[(^)'o^](r;o^) 

\E h \ < Ch*- 3 ^*- 1 + \D h \i) < c(h 2a - 4 + ^ 



< CK 



(3.38) 



where the latter term tends to zero owing to (3.26). Furthermore, we can prove that 
the second term in (3.35) is equal to zero. Indeed, let 
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It is easy to verify that tp h £ C/}(R) and tp h (0) = for every h > 0. Therefore, by 
(3.21) we obtain 

By (3.35) and (3.38) we conclude that 

J (R h E h (R h ) T )u 1 o 4) [($)' o y*](r, o -> 0. 

It remains to study the third integral in (3.32), which can be written as 



(3.39) 



/ (R h E h (R h ) T )i 
Jn 

-J 



fh rphf nh\T 



rd6 h 








\~dT 




\ V 


oy\ 


A3 


h 


^ V 


°y t l 


A3 


h 



dx 



/<W_ 
\dt 



°4)- 



dx 



We claim that 



dx = 0. 



(3.40) 
(3.41) 



To prove it, we consider again the sets Dh defined in (3.36). From (3.7), (3.31), 
(3.37), and from the boundcdncss of n we obtain 



/ 



(R h E h (R h ) T ) 13 



h 



"■ wd9 h 



— — O Zo 

dt 3 / 



(R h E h (R h ) T ) 1 



V dt 

v°yi 



l 



dx 



hu h J 



and the latter is infinitesimal owing to (3.26), so that (3.41) follows. In conclusion, 
combining (3.32), (3.34), and (3.39)-(3.41) we deduce that 



/i->0 



lim / (R n E n (R ,i y) 13 - 



[ E^rf 
Jn 



dx 



(3.42) 



for every 77 £ C 1 (R) with compact support and such that 77(0) = 0. Choosing a test 
function of the form 

^(x)=6 h ^-^(x 1 )) V (x 1 )e 1 , 
one can prove analogously that 



lim / {R h E h {R h ) T ) 12 ^±r 



h->0 



h 



E\\r( dx. 



(3.43) 



Let now Lp k £ C 2 (M.) with compact support be such that (fik(0) — ¥^(0) = 
for k = 2,3. We choose 77 = tp' 3 in (3.42) and 77 = tp' 2 in (3.43) and we add the 
two equations. Comparing with (3.25) and using the fact that E h (and therefore, 
R h E h (R h ) T ) is symmetric, we conclude that 

Entp'z + Eutp 3 ' + f 2 ip2 + hW3 dx\ = 



for every iff. £ C 2 (R) with compact support and such that (fk{0) = f' k (0) 
k = 2,3. By approximation we obtain (2.9) and (2.10) for all a > 2. 

Step 7. Euler-Lagrange equation for the twist junction 
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To conclude the proof of the theorem, it remains to verify the limiting Euler- 
Lagrange equation (2.11). We define 



where rj e C 1 (R) with compact support, 77(0) = 0, and 9 h is as in Step 6. Using <fr h 
as test function in the Euler-Lagrange equations (3.12), we obtain 



{R h E h (R h ) T ) 21 (6 h oz 3 l ) - (R h E h (R h ) T ) 31 (0 h oz%) (rj'oy^dx 
+ jf (R h E»(R>rhi o **) ((g)' o y?)^ dx 

l (R h E h (R h fhi o **) ((g)' o y^)^A dx 
+ jf [(R h E h (R h f) 32 (^oz%) {R*E\R h fh^o; 
+ h [ [f 2 (9 h o 4) - f 3 (8 h o 4)} ( V o y1) dx = 0. 



dx 



(3.44) 



Arguing as in the proof of (3.34), we can show that the first integral in (3.44) 
satisfies 



lim 

h-tO 



{R h E h {R h ) T ) 21 {6 h oz%) - (R h E h (R h ) T ) 31 (0 h oz' 2 l ) {q'oy^dx 

L 



=j (-E12 + E W W dxi. 
Jo 

The proof of (2.11) is concluded if we show that all other terms in (3.44) converge 
to zero, as h — >• 0. The last integral in (3.44) is infinitesimal owing to the estimate 

h f \f2(6 h ozZ)-M0 h O z%)]{r,oy*)dx < Ch(\\f 2 \\ L 4zZ\\ L 2 + \\f 3 \\ L 4z%\\ L2 ) 
Jn 

< Ch, 



which follows from (3.29) and (3.14). 
As for the term 

■dB h 



J Q [(R»E h (R h ) T ) 32 (^oz%) (R h E h (R h F) 23 (^-oz%) 



dx, 



we remark that by the symmetry of R h E h (R h ) T it can be written as 

L^***M[("°*)->] + [H" •*)]}*■ 

Arguing as in the proof of (3.41), we obtain that the above expression tends to zero, 
as h — > 0. 

It remains to prove that 

lim J \{R h E h {R h ) T ) kl o z f) [(#)' o yflfa o tf) dx = 
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for k,j 6 {2, 3}, k ^ j. To this aim, we fix k 
integral as the sum of two terms 



2, j = 3 and we write the previous 



^ - (R h E h (R h f) 21 (— o zl) [(#)' o ] („ o yl) dx 

nh 



1 r / ^7/9 \ i 

-(R*£ k (/i k f) 21 [(— o4') - l][(£f)'°rf](. I orf)<ix 

gg^MM l^),, (3 , 5) 



where < e < a — 2. Arguing as in the proof of (3.41), we obtain that the first 
term is infinitesimal. To study the second term, we notice that, if {ip h ) C C^(R) is 
a sequence of functions such that ip^(0) = and H^^HioofR) < C for all h > 0, then 
the map tp h (xi)ej can be used as test function in the Euler-Lagrange equations 
(3.12) for every h > 0, and we have 

(R h E h (R h ) T ) jl 



J 

Jn 



in 

If we now choose 



-[^ h )'oyh]dx = f Vf^oy^dx <CV\\f 3 \\ LHn) ^{). 

(3.46) 



V^i) := 



rj(s)ds, 



then by (2.28) we obtain 

l.h 



\\ip n \\ L °° < Ch a - 2 - e \\r]\\ L i < C for all h > 0, 

so that by (3.46) also the last term in (3.44) is infinitesimal as h — > 0. This concludes 
the proof of (2.11) and of the theorem. □ 

We conclude the section with a lemma, which provides us with a characterization 
of the limiting strain. This result is contained in the proof of [19, Theorems 4.3 
and 4.4]. We present here a concise proof for the reader's convenience. 

Lemma 3.1. Let all the assumptions of Theorem 2.6 be satisfied and let (R h ) be 
the sequence of rotations of Theorem 2.5. For every h > let G h : £1 — > M 3x3 be 
defined by 

(R h ) T V h y h - Id 



G 



d 2 (3 



and let G be the weak limit of (G h ) in L 2 (f2,M 3x3 ) (which exists, up to subse- 
quences, by (2.15)). Then, there exist g <E L 2 (0,L) and (3 e i 2 (i7,R 3 ), with zero 
average on S and dk(3 G L 2 {VL, R 3 ) for k = 2,3, such that, if we define 

M((3) := (x 2 A'e 2 + x 3 A'e 3 

we have 

'symM(/3) + («' + W v 2? + K) 2 ])ei ® ei if a = 3, 
sym M(f3) + u'e\ <g> e\ if a > 3, 

> sym M((3) + ge\ <g> e\ if 2 < a < 3, 

where u, Vk, and A are the functions introduced in Theorem 2.6. 

Proof. For every h > we consider the function ^ h : Q, — > R 3 defined by 

1 



symG 



(3.47) 



V*(aO : = Tzly h ( x ) - hx 2 R h { Xl )e 2 - hx 3 R h ( Xl )e 3 ] 



for every x G fi. By (2.17) we have that 



dkj — 1 Gek for every k = 2, 3. 



(3.48) 
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Therefore, if we define f3 h := -f h — 7 , where j h is the average of j h on S, we deduce 
by Poincare-Wirtinger inequality that (3 h is uniformly bounded in L 2 (S1,R 3 ). It 
follows that there exists f3 e L 2 (f2,R 3 ), with zero average on S, such that, up to 
subsequences, (3 h f3 in L 2 (£l, R 3 ). Furthermore, by (3.48) we have that dk(3 = Ge^ 
for all k = 2, 3. 

As for the first column of G, wc remark that by (1.1) we can write 
R h G h e l = hd l7 h + ^{x 2 {R h )'e 2 + x 3 (R h )'e 3 ) - J^R^I 

= hd ^ h + J^2^( Rh y e2 + ^(R h )'^) - j ^la-i 1 ^ 1 dx2dx z- 

(3.49) 

By (2.17) we have that R h G h e 1 -± Ge x weakly in L 2 (ft,M 3x3 ). Moreover, by (2.15) 
there exists a function g e L 2 ((0, L), R 3 ) such that 

/ ~ dx 2 dx 3 g weakly in L 2 ((0, L), M 3 ), 

while (2.24) yields 

j^x 2 (R h )'e 2 + x 3 (R h )'e 3 -± x 2 A'e 2 + x 3 A'e 3 weakly in L 2 ((0, L), M 3x3 ). 

Finally, by the weak convergence of ((3 h ) in L 2 (f2,R 3 ) we have that hd\(3 h — > in 
W _1 > 2 (Q,R 3 ); thus, passing to the limit in (3.49) we conclude that 

G = (x 2 A'e 2 + x 3 A'e 3 + g | 9 2 /3 | d 3 /?) . 

To obtain (3.47) it is now enough to define 

fi:=P + x 2 (g ■ e 2 )ei + x 3 (g ■ e 3 )ei, 

so that 



symG = sym^ 2 A'e 2 + a; 3 A'e 3 + (g-ei)ei 

This concludes the proof for 2 < a < 3. For a > 3 a characterization of 5 can be 
given. Indeed, one can observe that 

iy h - R h ei ^ f (d lV h - 1) + (1 - R^) 



f d lV h - R h e x f 



hal dx 2 dx 3 

= (u h )' - h a - 3 sym(R h - Id) u , 

where (u h ) is the sequence introduced in (1.3). By (2.20) and (2.25) we obtain the 
thesis for a > 3. □ 
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